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Abstract 

We study the crossover behaviors that can be observed in the high- 
temperature phase of three-dimensional dilute spin systems, using a field- 
theoretical approach. In particular, for randomly dilute Ising systems we 
consider the Gaussian-to-random and the pure-Ising-to-random crossover, de- 
termining the corresponding crossover functions for the magnetic susceptibil- 
ity and the correlation length. Moreover, for the physically interesting cases 
of dilute Ising, XY, and Heisenberg systems, we estimate several universal 
ratios of scaling-correction amplitudes entering the high-temperature Wegner 
expansion of the magnetic susceptibility, of the correlation length, and of the 
zero-momentum quartic couplings. 

PACS Numbers: 64.60.Ak, 75.10.Nr, 75.10.Hk 
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I. INTRODUCTION. 



The critical behavior of randomly dilute magnetic materials is of considerable theoretical 
and experimental interest [1-5] . A simple model describing these systems is provided by the 
Hamiltonian 

<ij> 

where the sum is extended over all nearest-neighbor sites, Sj are M-component spin variables, 
and Pi are uncorrelated quenched random variables, which are equal to one with probability 
p (the spin concentration) and zero with probability 1—p (the impurity concentration). For 
sufficiently low dilution 1—p, i.e. above the percolation threshold of the spins, the system 
described by the Hamiltonian Tip undergoes a second-order phase transition at Tc{p) < 
Tcip^l). 

The nature of the transition is rather well estabhshed. In the case of the random Ising 
model (RIM) corresponding to M = 1, the transition belongs to a new universality class 
which is distinct from the Ising universality class describing the critical behavior of the pure 
system. This has been clearly observed in experiments [3] on dilute uniaxial antiferromag- 
nets, such as Fej;Zni_2.F2 and Mna.Zn1_a.F2, in the absence of magnetic field [6] and in Monte 
Carlo simulations of the RIM, see, e.g., Refs. [7-10]. The critical exponents are independent 
of the impurity concentration and definitely different from those of the pure Ising univer- 
sality class. Field-theoretical (FT) studies [11-16] confirm these results. The fixed point 
(FP) related to the pure Ising universality class is unstable with respect to the addition of 
impurities and the renormalization-group (RG) fiow is driven towards a new stable random 
FP that controls the critical behavior. 

Unlike Ising systems, multicomponent 0(M)-symmetric spin systems do not change their 
asymptotic critical behavior in the presence of random impurities. Indeed, according to the 
Harris criterion [17], the addition of impurities to a system which undergoes a continuous 
transition does not change the critical behavior if the specific-heat critical exponent a of 
the pure system is negative, as is the case for any M > 2. From the point of view of RG 
theory, the Wilson- Fisher FP of the pure 0(M) theory is stable under random dilution. 
The presence of impurities affects only the approach to the critical regime, giving rise to 
scaling corrections behaving as \t\^'^, where r is the reduced temperature and Ai = —a. 
The exponent Ai is rather small for the physically relevant cases M = 2 and M = 3 — 
a = -0.0146(8) (Ref. [18]) and a = -0.1336(15) (Ref. [19]), respectively— giving rise to 
very slowly decaying scaling corrections. Experiments on ^He in porous materials [20,21] 
and on randomly dilute isotropic magnetic materials, see, e.g., Refs. [22-24], show that the 
critical exponents of XY and Heisenbcrg systems arc unchanged by disorder (see also the list 
of results reported in Ref. [4]). But, in order to observe the correct exponents in magnetic 
systems in which the reduced temperature is usually not smaller than 10"^, it is important 
to keep into account the scaling corrections in the analysis of the experimental data [22] . 

In this paper we study the crossover behaviors that can be observed in the high- 
temperature phase of three-dimensional dilute spin systems. First, we consider the crossover 
from the Gaussian FP to the stable FP of the model, i.e. the random FP for M = 1 and the 
pure 0(M)-symmetric FP for M > 2. Such a crossover can be observed at fixed impurity 
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concentration by varying the temperature. If \T — Tc\/Tc > G, where G is an appropri- 
ate Ginzburg number [25], fluctuations are irrelevant and mean-field behavior is expected, 
while for \T — T^/Tc < G the asymptotic critical behavior sets in. This crossover is not 
universal. Nonetheless, there are limiting situations in which the crossover functions become 
independent of the microscopic details of the statistical system: This is the case of the crit- 
ical crossover limit of systems with medium-range interactions, i.e. of systems in which the 
interaction scale is larger than the typical microscopic scale [26]. In this limit the crossover 
functions can be computed by using FT methods: for 0{M) models precise results have been 
obtained in Ref. [27] by using the three-dimensional massive scheme and in Refs. [28,29] by 
using the minimal-subtraction scheme without e expansion. 

In Ising systems there is also another interesting crossover associated with the RG fiow 
from the pure Ising FP to the random FP. When the concentration p is close to 1, by 
decreasing the temperature at fixed p, one first observes Ising critical behavior, then a 
crossover sets in, ending with the expected random critical behavior. In a suitable limit in 
which p — > 1 this crossover is universal. The corresponding universal crossover functions can 
be computed by using FT methods. 

These crossover behaviors are investigated here by using the fixed-dimension perturbative 
approach in powers of appropriate zero-momentum quartic couplings. We determine the RG 
trajectories and the crossover functions of the magnetic susceptibility x the second- 

moment correlation length ^, defined from the two-point function 



= (poPx So • Sx), (1-2) 

where the overline indicates the average over dilution and ( ) indicates the sample average at 
fixed disorder. This study allows us to compute the corresponding effective exponents and 
to determine several universal ratios of scaling-correction amplitudes entering their high- 
temperature Wegner expansions. Beside x ^-iid ^, we also consider zero-momentum quartic 
correlations and appropriate combinations that have a universal high-temperature critical 
limit, such as 



M + 2 r^0+ ^3^2 



G22 = - lim -pr-^, (1.3) 



where r is the reduced temperature, X4 is the zero-momentum four-point connected corre- 
lation function averaged over dilution, i.e., setting /i = J2xPxSx, 



Vx. = {{l^-l^f)-^{l^-l^)\ (1.4) 



and X22 is defined by 



Vx22 = {li-li?-{li-lA . (1.5) 
Their high-temperature Wegner expansion is given by 
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X = XrT-^ (l + Xr,l^^^ + Xr,2T^^ + ...) , (1.6) 
e = irT-" (l + er,lT^^ + er,2T^^ + ...), (1.7) 

G# = Gl(l + G#,r,ir^^ + G#,r,2T^' + ...) , (1.8) 

where Ai 2 are the exponents associated with the first two independent scahng corrections. 
For dilute Ising systems, a recent Monte Carlo study [8] provided the estimate Ai = 0.25(3); 
a rough estimate of A2 is A2 = 0.55(15), of. Sec. IIIB. For XY and Heisenberg systems 
Ai = —a, while A2 coincides with the leading correction-to-scaling exponent of the pure 
model, A2 = 0.53(1) for M = 2 and A2 = 0.56(2) for M = 3, cf. Ref. [4]. The ratios Cr,i/XT,i 
and XT,i/G^^r,i for i = 1,2 are universal. Their determination may be useful for the analysis 
of experimental or Monte Carlo data. In Eqs. (1.6-1.8) we only report the leading term for 
each correction-to-scaling exponent, but it should be noted that there are also corrections 
proportional to r^^% r^^^, etc., that may be more relevant — this is the case of systems with 
M > 2 — than those with exponent r^^ . 

The crossover behavior in dilute models was already studied in Refs. [30,14] in the Ising- 
like case and in Ref. [31] for multicomponent systems. However, Refs. [30,14,31] studied 
the crossover and computed the related effective exponents with respect to the RG flow 
parameter, while we compute effective exponents with respect to the reduced temperature, 
which have a direct physical interpretation. 

The paper is organized as follows. In Sec. 11 we discuss the FT approach. We first intro- 
duce the effective Landau- Ginzburg- Wilson 0^ Hamiltonian and some general definitions. 
Then, we generahze the approach of Ref. [27] by showing how to compute the crossover 
functions of the magnetic susceptibility and of the correlation length in terms of an effective 
temperature. These exact expressions allow us to determine the temperature dependence of 
several quantities near the critical point and, as a consequence, to compute some universal 
ratios of scaling-correction amplitudes entering the high-temperature Wegner expansion of 
X, ^, G4, and G22 for dilute Ising, XY, and Heisenberg systems. These results are pre- 
sented in Sec. HL Finally, in Sec. IV we extend the computation to the whole crossover 
regime, determining RG trajectories and effective exponents for Ising, XY, and Heisenberg 
systems with random dilution. In the case of Ising systems, we also discuss the Ising-to-RIM 
crossover, give analytic expressions for the crossover scaling functions — details are reported 
in App. B — and explicitly compute the crossover function associated with the magnetic sus- 
ceptibility. In App. A we prove some useful identities among the RG functions introduced 
in the FT approach. 



II. RG TRAJECTORIES AND CROSSOVER FUNCTIONS 

A. Definitions 

The FT approach is based on an effective Landau-Ginzburg- Wilson Hamiltonian that 
can be obtained by using the rephca method [32-35], i.e. 

Hmn -Jd^xlY,^ [{d,<Pair + + E ^ («o + VoS,,) ct>lct>l \ , (2.1) 

I ia ijab ' J 
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(a) ISING-LIKE 



(b) M>1 




FIG. 1. Sketch of the RG flow in the coupling plane {u,v) for (a) Ising (M = 1) and (b) 
M-component (M > 1) randomly dilute systems. 

where a,b = 1,...M and i,j = 1,...N. In the limit N ^ the Hamiltonian H-mn with 
Mo < and fo > is expected to describe the critical properties of dilute M-component 
spin systems. Thus, their critical behavior can be investigated by studying the RG flow of 
H-MN in the limit A/" — > 0. For generic values of M and N , the Hamiltonian TImn describes 
M coupled A/"- vector models and it is usually called MN model [1] . I-Lmn is bounded from 
below for Nuq + Vq> and -Uq + > 0. But, as discussed in Ref. [36], in the limit — > 
the only stability condition is Vq > 0. Figure 1 sketches the expected flow diagram in the 
quartic-coupling plane, for Ising (M = 1) and multicomponent (M > 2) systems in the limit 
— > 0. The relevant region for dilute systems corresponds to u < and thus the relevant 
stable FP is the random FP (RIM in Fig. 1) for M = 1 and the 0(M) FP for M > 2. 

The most precise FT results have been obtained in the framework of the fixed-dimension 
expansion in powers of zero-momentum quartic couplings. In this scheme the theory is renor- 
malized by introducing a set of zero-momentum conditions for the one-particle irreducible 
two-point and four-point correlation functions: 



(2.2) 



where 5, 



ai,bj 



j = Sab^ij, and 




Z. m{uSai ,bj,ck,dl + vCai,bj,ck,dl) , 



(2.3) 



where 



Sai,bj,ck,dl — g {Sai,bj^ck,dl ~^ ^ai,ck^bj,dl + ^ai,dl^bj,ck) ; 

Cai,bj,ck,dl — g Sij^ik^il {Sabred + ^ac^bd + Sad^bc) ■ (2.4) 

In addition one defines the function Zt through the relation 
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The critical behavior is determined by the stable FP of the theory, i.e. by the common 
zero u*, V* of the /3-f unctions 



(3u{u,v) = m— 



Uo,Vo 



(Jv{u,v) = m— 



(2.6) 



whose stability matrix has positive eigenvalues (actually a positive real part is sufficient). 
The critical exponents are obtained by evaluating the RG functions 



91nm ' 



Vt{u,v) 



dlnZt 
dlnm 



at u*,v*: 



(2.7) 



(2.8) 



The six-loop pertubative expansions of the (3 functions and of the critical exponents are 
reported in Refs. [11,37]. 

In the MN model, the RG functions satisfy a number of identities. Along the u = axis 
we have 





dPu 




dp. 




^ dp, 
du 






du 




dv 




u=0 




dVcp 
du 


u=0 


dri^p 
dv 


u=0 


= 0, 




axis we obtain 








dPu 




d(3y 




MN + 2 


dPu 


du 


v=0 


dv 




M + 2 


dv 


du 


v=0 


MN + 2 
M + 2 


dv<p 
dv 


= 0, 

v=0 


dr]t 
du 


v=0 


MN + 2 drjt 
M + 2 dv 


= 0. 

v=0 



= 0, 



(2.9) 
(2.10) 



= 0, 



v=Q 



(2.11) 
(2.12) 
(2.13) 



These identities can be proved order by order in the pertubative expansion, see App. A. 
The second set of relations was already reported in Ref. [38] for M — 1. 

In the limit ^ 0, the perturbative expansions in powers of u and v are not Borel 
summable at fixed ratio u/v (Ref. [39] shows it explicitly for the zero-dimensional theory 
with M = 1, but the argument has general validity), except when u = that corresponds to 
the 0(M)-symmetric 0^ theory. For M > 2, this is a minor problem since the relevant FP is 
the 0(M)-symmetric one. On the other hand, this is a notable limitation of the perturbative 
approach for the RIM. Nevertheless, rather reliable results for the critical exponents of 
the RIM universality class have been obtained from the analysis of properly resummed 
perturbative series. Several methods have been used: the Pade-Borel method at fixed u/v 
or the strictly related Chisholm- Borel method, the direct conformal-mapping method, an 
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expansion around the Ising FP [7], the double-Pade-Borel and the conformal-Pade-Borel 
method [11]. which, at least in zero dimensions [40], are able to treat correctly the non- 
Borel summability of the expansions at fixed u/v. The FT estimates of the critical exponents 
obtained from the analysis of the six-loop expansions reported in Refs. [37,11] depend only 
slightly on the resummation method. For instance, Ref. [11] reports v — 0.673(8) and r] — 
0.029(3) from the direct conformal-mapping method, and v = 0.678(10) and r] = 0.030(3) 
from an analysis that follows the ideas of Ref. [40]. A second source of uncertainty is the 
position of the FP. Monte Carlo [7] simulations give u* = —18.6(3) and v* = 43.3(2), which 
are significantly different from the FT estimates [11] u* — —13(2) and v* — 38.0(1.5), 
obtained from the numerical determination of the stable common zero of the /9-functions. 
However, as discussed in Ref. [7], the critical-exponent estimates show a relatively small 
dependence on the position of the FP. By using the Monte Carlo results for the location of 
the FP in the u-v plane, one obtains [7] u = 0.686(4) and r] = 0.026(3), which are close to 
the above-reported ones, obtained by using the field-theoretical estimates of the FP. In any 
case, it is reassuring that the FT results are in satisfactory agreement with the Monte Carlo 
estimates of the critical exponents, i.e. [7] u = 0.683(3) and rj = 0.035(2). The comparison 
of the different analyses shows that all different resummation methods give results of similar 
accuracy. In particular, the more sophisticated analyses suggested in Ref. [40] and employed 
in Ref. [11] apparently do not provide more accurate results than those at fixed u/v. For 
this reason, in the following we only use the Pade-Borel and the conformal-mapping method 
at fixed u/v. In the latter case, for the singularity of the Borel transform we use the naive 
analytic continuation for iV ^ of the result for the cubic model reported in Ref. [37]. 
The results of Ref. [39] suggest that this should allow us to take into account the leading 
divergent behavior of the series at least for sufficiently small \u/v\ (in zero dimensions for 
-1/2 < u/v < 0). 



B. Renormalizat ion-group trajectories 



The RG trajectories in the plane {u, v) are lines which start from the Gaussian FP located 
at u = V = and along which the quartic Hamiltonian parameters Uq and Vq are kept fixed. 
They are implicitly characterized by the equation 

uZu{u,v) Uq 



F{u,v) 



s. 



RG trajectories can also be determined by solving the differential equations 

-X^ = (3u(u(X),v(^)), 

(11) 

-A- = /3.(^A),^(A)), 

where A e [0, oo), with the initial conditions 

m(0) = w(0) = 0, 
du 
dX 



A=0 



dv 
dX 



1. 



(2.14) 



(2.15) 



(2.16) 



A=0 
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The solutions ^(A, ,s) and v{X, s) provide the RG trajectories in the {u, v) plane as a function 
of s. The RG trajectories relevant for dilute spin systems are those with s < 0. The 
attraction domain of the stable FP is given by the values of Uq and Vq corresponding to 
trajectories ending at the stable FP, i.e. trajectories for which 



u{X — oo, s) — u*, v{X — oo, s) — V*. 



(2.17) 



The crossover functions from the Gaussian to the Wilson-Fisher stable FP have been much 
studied in the case of the 0(M)-symmetric theories, both in the field theory [27-29] and 
in medium-range models [26]. In order to determine the crossover functions along the RG 
trajectories, and in particular those related with the correlation length ^, the magnetic 
susceptibility %, and the reduced temperature r oc r — Tc, we extend the method of Ref. [27] 
to Hamiltonians with many quartic parameters, such as Hmn- Using the relations 



p(l,2) 
ai,bj 



1 

X = J^X 

.(2) 



(0) 



dr 



Sni 



ai,bj 



dr 



Uo,Vo 



A 7-1 



(2.18) 



and Eq. (2.7), and defining 



fl<t>{\ s) = V<t^{u{X, s),v{X, s)), 

f/t(A, s) = rjt{u{X, s),v{X, s)) - r]^{u{X, s),v{X, s)), 



we derive the following expressions 

f(A, s) = ^vo = A, 
X(A, s) = xvl = A^ exp 



X 



t(A, s) = t/vI = / 

J A 



°° 2 - fj^jx, s) 
ax exp 

A X"^ 



dy 



fit{y,s) 



(2.19) 

(2.20) 
(2.21) 

(2.22) 



where ^, x, and f are dimensionless quantities. One can easily verify that in the Gaussian 
limit, i.e. for A — > 0+ or f — > oo, we have u,v — 0(A), f}^{X,s) = O(A^), fit{X,s) = 0(A); 
therefore f^^ — > 1 and — > 1, as expected. 

Eqs. (2.20-2.22) ahow us to compute ^ and x as functions of f and s. We can then define 
effective exponents by taking logarithmic derivatives of ^ and x fixed s: 



i^efr(T, S) 



d\n.T 



7eff(T,Sj = - 



d\n.x 



Vesir, s) = 2- 



dlnx 



(2.23) 



One can easily check that ?7cff = 2 — '~fcs/^cS = V<t>- ^he other hand, 7cff 7^ '~f{u,v) and 
^eff 7^ z^(m, f) where ^{Ujv) and u^UjV) are the RG functions associated with the exponents 
7 and u. 
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III. UNIVERSAL RATIOS OF SCALING-CORRECTION AMPLITUDES 



A. General results 



In order to determine the scaling-correction amplitudes, we compute the crossover func- 
tions close to the critical point, i.e., for A — * oo or f — 0"*". For this purpose, we consider 
the expansion of the RG functions around the stable FP {u*,v*). We write 

Pu{u, v) ^ buu{u - u*) + buv{v - V*), 

I3y{u, v) ^ byuiu - u*) + byy{v - V*). (3.1) 

Then, using Eq. (2.15) we have the following behavior, in the limit A ^ oo and for values 
of s in the attraction domain of the stable FP, 

«(A, S)^U*+ «A,l(5)A-'^^ + «A,2(S)A-'^^ + 

^;(A, s)^v* + vx,i{s)X-'^' + v^,2{s)X-^' + (3.2) 
where u!i,u!2 are the eigenvalues of the matrix 

b b r ^ ^ 

and we arc keeping only the leading terms in powers of \~^^ and A"*^^. In particular, we 
have neglected terms of order A"^'^^, A~^^% etc., which may be as important as those of 
order A"'^^. Moreover, we have 



D _ «A.i(=s) ^1 - b 

Ki = — 



vv 



v\,i{s) h,u ui-b 



uu 



D _ ^A,2(s) UJ2-byy buv 

VX,2{S) byu UJ2 - buu 

These ratios are independent of s, as expected because they are universal. Indeed, as we 
shall see, they can be related to the universal ratios G22,T,i/(j4,T,i of the scaling-correction 
amplitudes of G4 and G22, cf. Eqs. (1.3) and (1.8). 

We also expand the RG functions associated with the critical exponents, 

uiu, v) = — ] + Vu{u - u*) + Uy{v - V*), 

7(m, v) = [2- r]^{u, v)] z/(m, ti) Pa 7 + 7„(m - u*) + 7^(f - v*). (3.5) 
and define the s-independent quantities 

7A,i = 7uRi + Iv, ^X,i = ^uRi + ^v, \ = t^iZ/, (3.6) 

for i = 1, 2. Then, using Eq. (3.2), we find 
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and 



X(A, s) = Xa(s)A"-" fl + E Xa,,(s)A-"" + ..."), 

Xa.(^) = ^(^-^)m^), (3.7) 



t(A, = T,(5)A-^/^ (^1 + X: rA,.(s)A-'^' + ... j , 

(3.8) 

Using Eqs. (3.7) and (3.8), we can derive the Wegner expansion of ^, %, and of the zero- 
momentum quartic couphngs u and v in terms of the reduced temperature r. We obtain 

e(T, S) = ^ris)?"" (l + j^ er.(s)T^^ + ... j , 

Us) = TxisY, U{S) = lyTx,^is) Vx,^is) rA(s)-^% (3.9) 



X 



and 



X(r, s) = Xr{s)f-^ (^1 + i:Xr.(s)r^» + ...j , 



Xr(s) = Xx{s)Tx{sy, Xr,i{s) = ~ A- (/+ A-) ^^''^^^^^ ^^'^^^ 

and also 

2 

V^f,s)^V* + Y,Vr,i{s)r^' + ..., Vr,iis)^Vx,i{s)Tx{s)-^% (3.11) 
j=l 
2 

u{f,s)^U* + ^Ur,i{s)f^' + ..., Ur,i{s) ^ RiVx,i{s)Tx{s)~'^\ (3.12) 
i=l 

The results of Ref. [41] allow us to identify 

G4{f,s) = v{t,s), G22{f,s) = ^u{f,s), (3.13) 

and to obtain the corresponding scaling-correction amplitudes G4,r,i and G22,T,i defined in 
Eq. (1.8). 

From the above-reported relations we derive the following expressions for the universal 
ratios of scaling-correction amplitudes 
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— ^i, 



(3.14) 



X.T,i I'X.i T 



Their universality is explicitly verified since they are independent of s = uq/vq- 



B. Results for dilute Ising systems 

Using the results reported in Sec. Ill A, we can estimate several universal scaling- 
correction amplitude ratios. We analyze appropriate perturbative series that can be derived 
from those of the j3 functions and the critical exponents. Again, we use the conformal- 
mapping method and the Pade-Borel method at fixed ratio u/v. The errors we report take 
into account the resummation error and the uncertainty in the location of the FP. We com- 
pute each quantity at the FT and at the Monte Carlo FP. The final error is such to include 
both estimates. 

As a first step in the analysis we computed the subleading exponents and the ratios Ri 
and i?2- The exponent uji was already computed in Ref. [11], obtaining uji = 0.25(10) (using 
the double- Pade-Borel and the conformal- Pade-Borel method) and loi — 0.34(11) (using the 
direct conformal-mapping method), in substantial agreement with the Monte Carlo result 
u!i = 0.37(5) of Ref. [8]. In those analyses the field-theoretical estimates of the FP was 
used. We tried to compute uji by also using the Monte Carlo estimate of the FP. However, 
all methods gave largely fiuctuating results and no estimate could be obtained. Then, we 
determined UJ2- In this case, the conformal-mapping method provided reasonably stable 
results up to the Monte Carlo FP. We obtained [42] 002 = 0.8(2). 

Similar analyses were done for Ri and R2. Our final results are 

= -0.90(2), i?2 = -0.7(3). (3.15) 

Finally, we determined the ratios of scaling-correction amplitudes using relations (3.14). In 
order to have a check of the results, for each quantity we considered several series with the 
same FP value. We obtained 

Xr,i/&,i = 1.99(4), 
Xr,i/G4,.,i = -1.0(3), 

G22,T,l/Gi^r,l = 2.1(1), 
Xr,2/Cr,2 = 1.7(2), 

Xr,2/G4,.,2 = -0.4(2), 

G'22,r,2/G'4,r,2= 1.6(7). (3.16) 

The errors take into account the results obtained from different series and different resum- 
mation methods, and also the uncertainty on the location of the FP. It is interesting to 
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note that the rcsTilts for the ratios XT,i/^T,i show that the quantity x/^^ has much smaller 
scaling corrections than x ^"^^ ^- This fact was used in Ref. [7] in order to obtain a precise 
Monte Carlo estimate of r] from the high-temperature behavior of x/^^. For comparison, we 
report the corresponding values for the pure Ising universality class. Prom the analysis of 
high-temperature series one obtains Xr,i/Cr,i = 1.11(12) (Ref. [43]) and XT,i/Cr,i = 1.32(10) 
(Ref. [44]), while field theory gives [21] Xr,il^r,i = 1-47(4) and Xr,i/G'4,r,i = -0.30(4). 



C. Results for dilute multicomponent systems 

As in the Ising case, we determine the universal ratios of scaling-correction amplitudes 
by analyzing the six- loop expansions of the MN model [11]. Since the corresponding RG 
functions must be evaluated at the 0(M)-symmetric FP, i.e. along the u — Q axis, the series 
are Borel summable and the standard conformal-mapping method works well. 

Identity (2.9) allows us to obtain the following exact results for the universal quantities 

Ri- 

Ri = -1, R2 = 0, (3.17) 
which hold independently of M. We also obtain 

Xr,i/&,i = 1.97(2), 

Xr,i/G^,r,i = -17(2), (3.18) 

for dilute XY systems, and 

Xr,l/Cr,l = 1.97(2), 

Xr,i/G4,.,i = -2.5(4), (3.19) 

for dihite Heisenberg systems. The ratios Xr,2/^T,2 and XT,2/G4^r,2 are just the universal 
ratios of scaling-correction amplitudes of the 0(M)-symmetric models. Ref. [27] reports 
Xt,2/^t,2 = 1-57(2) and Xt,2/^t,2 = 1-63(4) respectively for XY and Heisenberg systems. 
We add the results XT,2/G^4,r,2 = —0.47(5) and XT,2/G^4,r,2 = —0.59(5) again for XY and 
Heisenberg systems. We finally mention an e-expansion study of the universal ratios of 
scaling-correction amplitudes [45], where the specific- heat and low-temperature quantities 
are considered. These results differ significantly from those determined in experiments [23] 
on Ni8o-xFe^(B,Si)2o- 



IV. CROSSOVERS IN RANDOMLY DILUTE SPIN SYSTEMS 

A. Crossover from Gaussian to random critical behavior in Ising systems 

In the case of the RIM, the FP's have been determined by using FT and Monte Carlo 
methods. For the random FP, we mention again the estimates u* = —18.6(3) and v* = 
43.3(2) obtained by Monte Carlo simulations [7] and the FT results reported in Ref. [11], 
u* — —13(2) and v* — 38.0(1.5). The position of the unstable Ising FP is ur — 0, vj — 
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23.56(2) (Ref. [43]). The RG trajectories for s > are not interesting for dilute systems; we 
only mention that they are attracted by another stable FP with O(A^) symmetry (A^ 0), 
located at [4,46] u = 26.3(4), v = 0. 

In Fig. 2 we show the RG trajectories for several values of s in the interval — 1 ^ s < 0, 
as obtained by numerically integrating the RG equations (2.15), after resumming the (3- 
functions. The figure has been obtained by using a single approximant, but others give 
qualitatively similar results. The resummation becomes less and less effective as |s| increases. 
This is expected since the singularities that make the perturbative series non-Borel summable 
play an increasingly important role as |s| gets larger. In any case, for — 1 ^ s < 0, the RG 
trajectories flow towards the random FP. For s ^ — 1, Pade-Borel resummations (in this case 
we cannot use the conformal-mapping method since the singularity we use is on the positive 
real axis) hint at runaway RG trajectories. If this is true and not simply an artifact of the 
perturbative approach, this suggests the existence of a value Smin ~ — 1 such that systems 
corresponding to s < Smin do not undergo a continuous transition. As a consequence, since 
u is directly related to the variance of disorder, the continuous transition is expected to 
disappear for sufficiently large disorder. This prediction may be checked by considering a 
lattice version of the continuum Hamiltonian 



(4.1) 



where ip is a scalar field and r{x) is a Gaussian uncorrelated random variable. Such a model 
is the starting point of the FT studies of dilute systems and, by using the replica trick, can 
be shown to be equivalent to the model with Hamiltonian (2.1). Our results suggest that 
there is a critical value Vc such that, if the variance of r{x) is larger than Vc, the continuous 
transition disappears. 
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FIG. 3. Ising systems: The effective exponents 7eff and v^s for several values of s in the interval 
-1< s < 0. 



Beside s = Smin, there is a second interesting value of s, the value s* such that the RIM 
FP is approached from above for s* < s < and from below for s^i^ < s < s*. One can 
easily realize that for this particular value of s the leading scaling corrections proportional 
to r^^ — and more generally proportional to r"'^^ — are not present in the Wegner expansions 
of the thermodynamic quantities. Numerically, by using the conformal-mapping method, 
we obtain s* = -0.25(5). 

The behavior of the RG trajectories close to the RIM FP can be determined by using 
the results presented in Sec. Ill A. We find that v{X, s) can be expanded as 

v{X,s)^v* + -{u-u*)-u,,,{s)(---) , (4.2) 

Ri \Ri R2J \ux,i{s)J 

where Ri and R2 are universal constants reported in Sec. Ill A, cf. Eq. (3.4), and ma,i(s) 
and ma,2("5) are expansion coefficients defined in Eq. (3.2). Note the presence of the non- 
analytic correction which shows that, close to the FP, trajectories are only defined for 
{u — u*)/ma,i(s) > 0. This is expected on the basis of general arguments [47-49]: along any 
RG trajectory one expects nonanalytic corrections proportional, for instance, to nui/uJi + m, 
n, m being integers. 

Using Eqs. (2.20), (2.21), and (2.22), one can compute the crossover functions ^ and 
X along the RG trajectories, i.e. for fixed s, and the corresponding effective exponents UeS 
and 7eff, cf. Eq. (2.23). The effective exponents 7eff and are shown in Fig. 3 for several 
values of s < within the attraction domain of the RIM FP. We note that they become 
nonmonotonic for s < s* ^ —0.25, where the RG trajectories reach the RIM FP from below, 
see Fig. 2. 
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The crossover from the Gaussian FP to the RIM FP has also been investigated in 
Rcfs. [30.14] in the framework of the minimal-subtraction scheme without e expansion. The 
effective exponents computed here differ from those defined in Refs. [30,14], since there the 
nontrivial relation between temperature and RG flow parameter was neglected. In spite 
of the different definitions, the crossover curves obtained in Refs. [30,14] have the same 
qualitative features of those shown in Fig. 3. 

The above field-theoretical results may be related with those obtained in a specific (lattice 
or experimental) system by comparing the behavior in a neighborhood of the critical point. 
Given a quantity O, we can write for the field-theoretical model 

{O) ^ C{s)t-'' (l + A(s)t^i + A2(s)t^^) , (4.3) 

while for the lattice or experimental system we write 

{O) « Dt-'^ (l + B,T^^ + B,T^-) . (4.4) 

Then, we require these two expansions to agree apart from a rescaling of the reduced tem- 
peratures Tg — cr, i.e. 

B, = c^Mi(s), B2 = c^'A^is), (4.5) 

which gives 




A1/A2 



(4.6) 



Thus, in order to match the two expansions one should first determine s by using Eq. (4.6) 
and then fix c by using Eq. (4.5). This provides a mapping between the field-theoretical 
model and the considered system. This relation does not depend on the chosen quantity O 
because of the universahty of ratios of subleading corrections (the ratios of the Si's and of 
the B2^s of two different quantities are universal). Note, however, that the existence of this 
mapping is not guaranteed. In particular, Eq. (4.6) requires Ai{s)/Bi and A2{s)/B2 to be 
both positive. Since Ai{s) changes sign for s = s*, it is always possible to have Ai{s)/Bi > 0. 
But there is no guarantee that A2{s) / B2 can always be made positive. This is the well-known 
sign problem that has been discussed at length in 0{M) models [50-52,48]. For instance, 
it prevents to match the crossover curves for the scalar 0^ theory with the results obtained 
for the three-dimensional Ising model. Ref. [52] suggested the use of the "strong-coupling" 
branch g > g*, but this proposal fails in the massive zero- momentum renormalization scheme 
because of the nonanalyticity of the RG functions at the FP [47-49] . This phenomenon is 
even more evident in the RIM case, cf. Eq. (4.2). It should also be stressed that the 
mapping defined by Eqs. (4.5) and (4.6) does not imply that the field-theoretical crossover 
curves exactly match the corresponding ones for the considered system for all values of r. In 
particular, there is no relation among the neglected coefficients in the Wegner expansions. 

Finally, let us discuss the RIM with nearest-neighbor interactions on a cubic lattice 
with spin density p. Numerical simulations show the presence of a dilution-independent 
continuous transition up to p = 0.40 [8] . It is usually conjectured that the transition persists 
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up to the percolation threshold of the spins p = pc, Pc = 0.3116081(13) on a cubic lattice 
[53]. Below the percolation point the spins form finite domains and are therefore unable 
to show a critical behavior. It should be remarked that the transition for p = p^ is not 
described by the field-theory model (2.1) and thus the RIM for p = Pc does not correspond 
to s = Smin- For the same reasons, the fact that the transition disappears for p < Pc does 
not provide evidence in favor of a finite Smin- However, if the RIM can be related with 
the field-theory model (for example if Eq. (4.6) can be solved for any value of p) and the 
RIM with p ^ Pc corresponds to the field-theory model with s — > s, then we can conclude 
|s| < |>Smin|. Now we show that this condition is approximately verified. For this purpose, 
we must determine the relation between the RIM and the field-theory model. We use the 
results of Ref. [34] that map the RIM onto a translationally-invariant effective Hamiltonian 
-^RiM ^ field (f). The expansion of 7{^^^ for — has the same form, up to order 0'*, 
of the Hamiltonian (2.1) with M = 1. The corresponding quartic couplings u^^^ and v^^^ 
appearing in this expansion are related to the magnetic concentration p (note that such 
result does not depend on the lattice type and on the spin-spin interaction as long as it is 
of short-range type) by 

C™ocp(p-l), ^o^^^ocp, (4.7) 

and in particular 



u^'"^ 3 



,RIM 



2 



1-p). (4.8) 



It is tempting to assume s ~ u^^^'^/v^™, which means that we neglect the fact that in Ti^™ 
there are interactions 0" with any n > 4. The relation s ^ —3(1 — p)/2 follows. Using this 
relation and the numerical results of Refs. [8,7], we can get an independent approximate 
estimate of s*. Since in the RIM on a cubic lattice one does not observe the leading scaling 
correction for p* ^ 0.8, we obtain s* ~ —0.3, which is reasonably close to the FT estimate 
s* — —0.25(5). Moreover, the percolation threshold Pc — Pc = 0.3116081(13) on a cubic 
lattice [53] — apparently corresponds to s ~ —1, which is compatible with the predicted 
inequality |s| < 



'mm 



B. Crossover from Ising to random critical behavior 

The FT approach presented in Sec. 11 allows us to determine also the Ising-to-RlM 
crossover functions. Considering in general a quantity O that behaves at the Ising FP, 
i.e. in the absence of disorder, as t~'^'^' , standard RG arguments show that, in the limit 
p 1 and t = {T — Ti)/Ti — > 0, where Tj is the critical temperature of the pure Ising 
model, O can be written in the scaling form 

O - aot-P'^'Bo {gt-^) - a.e'Co {gC'^^"') , (4.9) 

where g cc 1—p is the scaling field associated with disorder, which is a relevant perturbation 
of the Ising FP, and ao and ai are normalization constants. The crossover exponent is 
equal [1] to the Ising specific-heat exponent aj, (f) — aj. The functions Bq and Co are 
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FIG. 4. Ising systems: The quantity x(A, s)A ^^^'^ as a function of \s\^°''^'^' for several values 
of s. 

universal, apart from trivial normalizations. By properly choosing qq and ai we can require 
-80(0) = Co{0) — 1- Another condition can be added by properly fixing the normalization 
oig. 

Within the FT approach the limit g ^ corresponds to s — 0~ and g^'^^'^' ~ sX'^^'^' . 
Therefore, crossover functions are obtained by taking the limit s — > 0~ and ^ = A — > 00 
of the quantity 0^~^', keeping sX°''^'^' fixed. In Fig. 4 we show numerically that such a 
limit exists for the susceptibility x- We consider x^'"^^^' = xi-^^ 3)X~'^^^' and plot this 
combination as a function of \s\^'^'^'^' for several values of s. The curves, obtained by using 
Eq. (2.21) and the conformal-mapping method, rapidly converge to a limiting function. 

In order to compute the crossover functions, we must first study the limit s ^ 0~ of the 
RG trajectories. As it can be seen from Fig. 2, in this limit the trajectory will eventually be 
formed by two parts connecting at the Ising FP: the line u = starting at the Gaussian FP 
and ending at the Ising FP, and a line v = g{u) connecting the Ising FP to the RIM FP. The 
line V = g{u) corresponds to a RG trajectory and therefore {u{X),v{X)) = [u{X), g{u{X))] 
must satisfy Eq. (2.15). Therefore, g{u) is the solution of the differential equation 



dg_ 
du 



(3y{u,g{u)) 
Pu{u,g{u)) 



(4.10) 



with the initial condition g(0) — vj. As discussed in App. B, g{u) is expected to be analytic 
for 14 —> and thus it can be expanded as 



g{u) =vi + J29n 



u 



(4.11) 
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FIG. 5. Plot of the curve v = g{u). The full line represents the quadratic curve given in 
Eq. (4.13), while the dotted line shows the linear approximation g{u) = vj — u. We also show the 
position of the RIM FP as obtained by MC simulations of the RIM (circle) and by FT calculations 
(triangle) . 

In App. B3 we compute the first coefficients: gi = —1, a consequence of identity (2.9), 



g2 = 0.0033(1), and gs = 1(2) x 10"^ 

Since g{u) corresponds to an RG trajectory with s = 0, Eq. (4.2) implies that, close to 
the RIM FP, we have 



Eq. (4.12) shows that g{u) is not analytic at the RIM FP. Of course, one should check that 
^A,2(0) does not vanish. We are not able to verify numerically this condition, but we believe 
that it is unlikely that ua,2(0) = 0. Indeed, the curve g{u) is a special curve only at the 
Ising FP, but it has no special status at the RIM FP and thus it should be nonanalytic as 
any generic RG trajectory [54]. 

The curve g{u) can be computed [55] by resumming the perturbative series for the f3 
functions and then by explicitly solving Eq. (4.10) with the initial condition g{0) = vj. The 
result turns out to be very well approximated by the simple expression 



where vj = 23.56(2) is the coordinate of the Ising FP [43] and g2 ~ 0.0033. Such an approx- 
imation is effective, within the resummation errors, up to the RIM FP. A graph is reported 
in Fig. 5. The results obtained by using [3,1], [4,1] and [5,1] Pade-Borel approximants would 
not be distinguishable from the curve (4.13) shown in Fig. 5. For instance, g{—13) ~ 37.1 



(4.12) 



g{u) ^ vj-u + g2U^, 



(4.13) 
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and — 18.6) ~ 43.3, so that Eq. (4.13) is perfectly compatible with the MC estimate of the 
FP, u* = -18.6(3), V* = 43.3(2), and with the FT result, u* = -13(2), v* = 38.0(1.5), see 
Fig. 5. The fact that both estimates lie on the limiting curve v = g{u) shows that the FT 
approach is effective in determining the Ising-to-RIM trajectory, although it is apparently 
unable to determine precisely the position of the FP on this curve. As a final check, we com- 
pute g'{u*). Using Eq. (4.12) and the estimate of Ri reported in Sec. IIIB, Ri = —0.90(2), 
we obtain g'{u*) = -1.11(2), while Eq. (4.13) gives g'{u*) = -1.12 (resp. -1.09) at the 
Monte Carlo estimate (resp. field-theoretical) of the FP. The agreement is satisfactory. 

Once we have determined g{u), we can compute u{X, s) in the crossover limit. In App. B 1, 
we show that, in the crossover hmit, s — > 0~ keeping \s\X°'^^'^' fixed, u{X,s) converges to 
U{a) which is implicitly defined by 



a — U [a) exp < / 

[ vi Jo 



dx 



1 ^ 



where 




(3u{x,g{x)) ajx 



(4.14) 



(4.15) 



and El and E2 are normalization constants such that t/(o") ~ a for a — > 0. Their explicit 
expressions are reported in App. B 1. Of course, v(A, s) = giJJioy) in the scaling limit 
s — > 0~. The curve g{u) and Eq. (4.14) completely fix the relevant RG trajectory in the 
crossover limit. 

The computation of the crossover functions is then completely straightforward. We 
consider the RG function 0(A, s) associated with O and assume that it satisfies the RG 
equation 



( \n 



(4.16) 



where p{u,v) is the corresponding RG function such that p{0,vi) = p/. The crossover limit 
is studied in detail in App. B2. We find that the crossover function Coiy) can be written 
as 



Co{y) = exp 



ru{<y) p{x,g{x))-pi 
Jo I3u{x,g{x)) 



(4.17) 



where the relation between y and a should be fixed by choosing an additional normalization 
condition. 

We wish now to specialize the previous discussion to the magnetic susceptibility. In this 
case p{u,v) = 2 — ri^{u,v). In order to completely specify the function C^{y) appearing in 
Eq. (4.9) we must fix the normahzation of g. We use the small-y expansion of C^{y). Since 
C^(0) = 0, see App. B 2, we require 



y' 



n=3 



(4.18) 
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FIG. 6. The crossover function C^{y) normalized according to Eq. (4.18). The dashed hne 
represents the asymptotic behavior (4.20). The inset shows the small-y behavior: the expansion to 
order corresponds to the dotted hne while the expansion to order y'^ corresponds to the dashed 
line. 

for y ^ and C^{y) to be defined for y > 0. With these normalizations we have 



where y = —yocr. The constant yo is positive and is computed numerically in App. B3: 
yo = 0.072(8). The scaling function C^{y) is shown in Fig. 6. 

We study the small-y and large-y behavior of C^{y). A rough estimate of the coefficient 
C3 is C3 = —4(2), see App. B3. For large values of y, we have 



where rj is the RIM exponent. The best estimates of the exponents rjj and rj of the Ising 
and RIM universality classes are r]j = 0.03639(15) (Ref. [43]), 7]i = 0.0368(2) (Ref. [56]), 
and Tj = 0.035(2) (Ref. [7]). These results suggest rji > rj. This is confirmed by the analysis 
of the fixed- dimension FT series: all analyses find rji > rj. In particular, analyses based 
on an expansion around the Ising FP [7] find rji — rj = 0.002(2). This suggests that C^{y) 
diverges for large y with a very small exponent, 1^7(77/ — T])/aj = 0.01(1). We also estimated 
the coefficient Coo appearing in the large- ?/ behavior of C^{y), obtaining Coo = 1.05(5). We 
proceeded as follows. First, for given approximants of the RG functions, we computed the 
exponents rjj, uj, and r], and the function C^{y). Then, we calculated C^{y)y~'^'^''T~^^^°'' and 
determined the constant Cqo from its large-y behavior. This procedure gave an estimate of 
Coo for a given set of approximants. The final result was obtained as usual, by comparing 
the results of different approximants and of series of different order. 
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FIG. 7. RG trajectories in the dilute XY model for — 1 < s < 0. 

C. Crossover in randomly dilute multicomponent spin systems 

In the case of multicomponent systems, the stable FP is the 0(M)-symmetric FP 
(0,w*). Precise estimates of v* have been obtained by employing FT and lattice techniques 
[4,18,19,46,57]: v* = 21.16(5) (FT) and v* = 21.14(6) (lattice) for the XY universality class, 
V* = 19.06(5) (FT) and v* = 19.13(10) (lattice) for the Heisenberg universality class. 

In Figs. 7 and 8 we show, respectively for XY and Heisenberg systems, the RG tra- 
jectories in the u, v plane for several values of s in the range — 1 < s < 0. Figs. 9 and 10 
report the corresponding effective exponents r/gfr and respectively for XY and Heisenberg 
systems. They are nonmonotonic. In particular, for s close to —1, rjes becomes negative for 
intermediate values of f. As in the Ising case, the resummations become less reliables — and 
again hint at runaway trajectories — for s ^ — 1. 

Finally, we mention that the RG trajectories and the effective crossover exponents of 
dilute Heisenberg systems have been recently investigated in Ref. [31], using a two-loop ap- 
proximation within the minimal-subtraction scheme without e expansion and neglecting the 
nontrivial relation between temperature and RG flow parameter. In spite of all simplifying 
assumptions, the results are in qualitative agreement with ours. Moreover, Ref. [31] dis- 
cusses crossover phenomena observed in experiments on isotropic magnets, showing several 
results for the effective exponents that are in qualitative agreement with the curves shown 
in Fig. 10. 
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FIG. 8. RG trajectories in the dilute Heisenberg model for — 1 < s < 0. 



APPENDIX A: SOME RELATIONS AMONG THE RG FUNCTIONS 



In this Section we prove identities (2.9) and (2.10) holding along the m = axis, and 
(2.11), (2.12), and (2.13) holding along the w = axis. 

Let us first prove the identities along the u = axis in the case M = 1; the extension 
to other values of M is straightforward. We consider a generic theory with fields and 
Hamiltonian density 

^ = \ Ud,<p^y + ru^p^^r + u^p^^r + ^ e c^bcd<pu'<p''<p^. (m) 

tj.,A ^ A ^- A ^- ABCD 

For Cabcd = the model is simply a collection of decoupled Ising 0"^ theories. In order 
to compute the corrections to first order in Cabcd, we consider the one-particle irreducible 
correlation functions of the fields expressed in terms of the bare couplings g and Cabcd 
and of the inverse susceptibility as effective mass (the results also hold for the massless 
theory in dimensional regularization) 

rA„...,A„ = (<^^\...,<^^"). (A2) 

Then, we prove that, if all indices are equal, 

rA,A,-,A = fig) + Caaaa^^ + OiC^) = fig + Caaaa) + OiC^). (A3) 

Using this relation, one can derive identities (2.9) and (2.10). Indeed, Eq. (A3) implies that 
(setting uq = uo/m and vq = vo/m) 
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FIG. 9. The effective exponents rjeg and z/gfj of the dilute XY model for —1 < s < 0. In 
the Gaussian limit = 1/2 and rj = 0, while in the Wilson-Fisher limit v = 0.67155(27) and 
rj = 0.0380(4) (Ref. [18]). 



Z^ = U{uo + vo) + 0{v^,), (A4) 
u + v = fu+v{uo + vo) + 0{v^), (A5) 
Pu + /3v = Muo + vo) + Oivl). (A6) 

To prove Eq. (A3), consider a generic diagram D contributing to the correlation function. If 
is used as effective mass or the mass vanishes and dimensional regularization is used, the 
diagram has the following properties: it does not contain tadpole subgraphs; given a vertex 
V, the subdiagram D/V obtained by deleting the lines going out of V may be disconnected, 
but each piece contains at least one external line. The contribution of the diagram D is 
the product of three factors: the first is the integral over the internal momenta, the second 
the symmetry factor, and the third one — we call it I{D)a,a,...,a — takes into account the 
interaction vertex 

Vabcd = —g^ABCD — Cabcd- (A7) 
Clearly, we are only interested in the last term which can be written in the form 

J^{D)ji^A,...,A = H^)a,a,...,a\c=o ^ X! H^/^)a,a,...,A;I,j,k,l\c=oCijkl 

veD 

= (-g)--ni-gr-'CAAAA, (A8) 

where n is the number of vertices of D. In the last step, we have used the two properties we 
have mentioned above: they guarantee that I{D/V)a,a,...,A;LJ,k,l = {—g)"'~^^AijKL, since 
for C = a connected diagram does not vanish only if the indices on the external legs are 
all equal. Eq. (A8) gives immediately Eq. (A3). 
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FIG. 10. The effective exponents ryeff and feg of the dilute Heisenberg model for —1 < s < 0. 
In the Gaussian limit v = 1/2 and r] = 0, while in the Wilson-Fisher limit v = 0.7112(5) and 
r] = 0.0375(5) (Ref. [19]). 



Identities (2.11), (2.12), and (2.13) along the v = axis can be proved in a similar 
fashion. Let us again restrict ourselves to the case M = 1, the extension to generic values 
of M being straightforward. Consider the Hamiltonian density 

n = I Ud,<P^r + IrU^^r + ^ Y^i^^my + i E (A9) 

^ fi,A ^ A ^- AB ^- ABCD 

where Cabcd is symmetric in all indices. For Cabcd = the model is simply an iV- vector 
0^ theory, where N is the dimension of the field. In order to compute the corrections to first 
order in Cabcd, we consider here a different set of correlation functions: 0(A^)-invariant 
(therefore there are no external indices) one-particle irreducible correlation functions of the 
fields and of any 0(A^)-invariant operator. Consider again a diagram D, a vertex V, and 
the interaction contribution I{D/V)ij^k,l for C = 0. Because of the 0{N) invariance, its 
symmetrized part is given by 

I{D/V) {I,J,K,L}\ ^ = I{D/V){6uSkl + SikSjl + SilSjk) ■ (AlO) 

Then, repeating the argument leading to Eq. (A8), we obtain 

m = imc=,-^T.hD/V)Y.Cnjj. (All) 

V IJ 

The constant Ylv 1^) is determined by computing the derivative of I{D) with respect 
to (? at C = 0. 
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= E E i{D/v) 

I,J,K,l\c=q ^ ( ~q ) i^IJ^KL + SikSjl + SilSjk) 
C=0 VeDIJKL ^ "^^ 

= - 5^ i{D/V)N{N + 2) . (A12) 
It follows 

/p) = /(^ + ^|ff^)+0(C^), (A13) 

where f{g) — 7(D) |^^q. This relation is valid only for 0(A'")-invariant quantities, but it can 
also be applied to the correlation functions of the elementary fields by simply contracting 
the external indices. It allows us to derive a number of relations involving the /5-functions 
and the RG functions associated with the exponents. For example, considering the MN 
model (2.1) for M = 1, relation (A13) implies Eqs. (2.11), (2.12), and (2.13) with M = 1. 

APPENDIX B: THE ISING-TO-RIM CROSSOVER 

In this appendix we compute the limit s ^ 0~ of the RG trajectories and the Ising-to- 
RIM crossover function Co{y), cf. Eq. (4.9). 



dI{D) 
dg 



1. The limit s ^ of the RG trajectories 

Here, we wish to prove Eqs. (4.14) and (4.15) that give u{X,s) in the crossover limit 
s — > 0~ keeping \s\X°'^^'^' fixed. As discussed in Sec. IV B, in the crossover hmit the RG 
trajectory is formed by two parts connecting at the Ising FP: the hue u — starting at the 
Gaussian FP and ending at the Ising FP, and the line v = g{u) connecting the Ising FP to 
the RIM FP. Now, we will solve the flow equations (2.15) in the two cases and we will match 
the two solutions in the neighborhood of the Ising FP. Let us consider first the behavior 
near v — g{u). The flow equation for u{\, s) can be written as 



-A^ = Pu{u,g{u)). 
Since Pu{u,vi) = —uai/vi for m — > 0, we can write the solution as 



(Bl) 



A = A{s)u{\ sf'"' exp 



u{\,s) 



dx 



+ 



^/3u{x,g{x)) ajx 



(B2) 



where A{s) is a (at this stage unknown) function of s. 

Now let us consider the second case, i.e. the trajectory near the u = axis. For u ^ 0, 
we can write Pu{u,v) = uf{v) + 0{v?), with /(O) = —1, /(f/) = —aj/ui. As for j3v{u,v) we 
simply set u = 0. Note that /3^(0, v) = —v + 0{v'^) for — >• and /5^(0, v) = —u!i{yi — v) for 
V ^ vj. In the hmit we are interested in, the RG equations (2.15) become 



-\— = P,{0,v), 



-A^ = ./(.). 



(B3) 
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Keeping into account the initial conditions (2.16), we obtain 



X — v exp 



— dx 



1 1 

+ 



r, ( fix) i\ 

u — sv exp / dx -— — — r . 

Jo \lJy[0,x) xj 



(B4) 
(B5) 



Eqs. (B4) and (B5) implicitly define u{X, s). We must now match the two solutions near the 
Ising FP, determining the unknown constant A{s). If we define 



El = vi exp / dx 
Jo 

T,2 = vi exp 



fix) 



^/3i,(0,,T) X UjUJiivj-x)^ 

1 



dx 



1 1 
+ - + 



^Pvi^^x) X ujiivi — x) ^ 
for V ^ vj Eqs. (B4) and (B5) can be written as 

Therefore, for i> — > we have 



(B6) 



ii(A,s) ft! sSi 



(B7) 



(B8) 



On the other hand, Eq. (B2) gives for m ^ 0, 

ii(A,s) ft (A/A(s))"^/'^^ (B9) 
By comparing Eqs. (B8) and (B9) we obtain Ais). Finally, Eq. (B2) can be written as 



ail VI 



m(A, s) exp 



ai 



«(A,s) 



dx 



+ 



Mx,gix)) ajx 



(BIO) 



This ends the proof of Eqs. (4.14) and (4.15). 



2. Crossover functions 

The computation of the crossover function is similar to that presented in App. B 1. We 
first consider the RG equation (4.16) on the line v — giu). Using the flow equation for 
ti(A, s) we can write 

do ^ _iM^a (BID 

du Puiu,g{u)) 

The solution can be written as 



26 



O = B{s)u{X, sY'"'/"' exp 



Jo \Pu{x,g{x)) ajx 



(B12) 



where B{s) is an unknown function. 

For — > 0, we can use the flow equation for v{X,s) and write 



dO 
dv 



■O. 



(B13) 



Wc assume p(0, 0) = po (po is the naive Gaussian dimension of O) and (9 ~ OqX^'^ at the 
Gaussian FP (Co is a normahzation constant). Then, the previous equation gives 



O = OqV^'" exp 



— dx 



p{0,x) ^ Po' 



10 \(3v{0,x) X 
Now, we must compute the behavior for v ^ vi. Defining 



(B14) 



Ti = Oo-yf exp 
we obtain for v ^ vj 



dx 



p(0,x) 



Pi 



(3y{0,x) X u!i{vi — x)^ 



(B15) 



vipi/ai 



(B16) 



where we have used Eq. (B7). On the other hand, Eq. (B12) gives in the hmit ti — > 



O ^ B{s)u'"P'l'^'. 



Therefore, 



= Ti 



u 



uipl/ai 



exp 



dx 



+ 



^(3u{x,g{x)) ajx 



(B17) 



(B18) 



FinaUy, by using Eq. (BIO) to ehminate u^ipiI°'i ^ we obtain 



O^TiEa^'A^'^exp 



L 



^^"^x P^^'3{x))-Pi 



The crossover function Co{y) normahzed so that Co(0) = 1 is then given by 



(B19) 



Co{y) = exp 



^^"^ p{x,g{x)) - pi 
I3u{x,g{x)) 



(B20) 



To fully specify the function Co{y) we must also relate y with a by adding an additional 
normalization condition. For the magnetic susceptibility this is done in detail in Sec. IV B. 

We can specialize these results to the observables we have considered in the paper. First, 
we consider the four-point quartic couplings G22 and G4. Since they are related to u and v. 
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= u/2) and G4 — v (see Ref. [41]), and u !v U{a), v g{u) in the crossover limit, we 
obtain 

CgM = ^, (B21) 

CgM = -9{U{a)). (B22) 

Vi 

Note that 6*022(2/) is not simply ^'"(cr) since the crossover function is defined by ~ 
\°"/^'CG22{y)- These equations can also be derived from Eq. (B20) bu using p{u,v) — 
—Pu{'u,v)/u and p{u,v) = —Py{u,v)/v for u and v respectively. 

Finally, let us consider the magnetic susceptibility x- In this case p{u,v) = 2 — rj^lujv). 
Thus, by using Eq. (B20) we obtain Eq. (4.19). Let us now show that 6*^(0) = 0. First, 
note that, because of identity (2.10), near the Ising FP we have 

r]^{u, v)-r]i = A{u + v- vi)[l + 0{u) + 0{v - vi)], (B23) 

where A is a constant. Then, since g{u) — vj — u + O^u^), we obtain 77<^(w, g{u)) —r]i — 0{u^). 
Substituting in Eq. (4.19), this gives immediately C*^(0) = 0. 

Finally, we argue that the crossover function Co{y) and g{u) (that can be related to the 
crossover function oi v — G4) are analytic for y — > and — > respectively. This is not 
obvious since for it = RG functions are nonanalytic at the Ising FP [47,49]. We will now 
show that such a problem does not arise for the RG functions defined along the crossover 
line V = g{u). The reason is that such a line has a very special status at the Ising FP: 
It is the line that is tangent to the relevant direction associated with disorder and that is 
orthogonal to all irrelevant directions. 

To clarify the issue, let us for instance consider the singular part of the free energy. In 
a neighborhood of the Ising FP it can be written as [58] 

Fsing = ft^'^F {f,f,-^, {fj,^^}) , (B24) 

where ft, fp, and {/j} are the nonlinear scaling fields associated with the temperature, the 
dilution, and the irrelevant RG operators. For t = {T — Ti)/Ti — >■ and p ^ 1, ft ^ t and 
fp ^ {1 — p) ^ g- The exponents Aj are associated with the irrelevant operators and are 
positive. A basic result of RG theory is that the nonlinear scaling fields are analytic in t and 
p and the function F is analytic in all its arguments. In the crossover limit, /j approaches 
a constant and ft goes to zero, so that fif^' 0. It follows 

i^sing~i'/'^^i^(^r^{0}), (B25) 

which shows that the crossover function associated with Fging is analytic in gt~'l'. The 
argument can be trivially generalized to any zero-momentum quantity; we conjecture that 
it also applies to quantities involving the correlation length. 

3. Some numerical results 

In this Section we report some details on the numerical computation of g{u) and C^{y). 
Let us first focus on the determination of the coefficients gn defined in Eq. (4.11). They 
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have been obtained by rcsumming pcrturbative series Qniv) such that Qn — gn{vi)- For the 
purpose of determining gn{v), we write 

f3uiu,v) = Y.b^Mu", (B26) 

n 

(3,{u,v) = Y,K,n{v)u\ (B27) 

n 

Then, by using Eq. (4.10), we obtain 

b'Uv) - - 2b^,2{v) 



92 = 



2[ai/ui + b^,i{v)] 



(B28) 



V=Vl 



and similar, but more complex, expressions for g^, g^, etc. The series gn{v) can be obtained 
by expanding the right-hand side in powers of v. For g2 and g^ we obtain 

g^(v) = 0.00663146 - 0.00693165^; + 0.0116887^^ - 0.0225971# + (B29) 

+0.0455962^;^ - 0.0954011t;^ + O(w^), 
g^{v) = 0.0000293176 - 0.0000813454^ + 0.000206937^;2 - 0.000485549^^ + (B30) 

+0.00110105v^ + O(v^), 

where v = 3i'/(167r). By resumming these series we get 

g2 = 0.0033(1), gs = 1(2) x 10"^ (B31) 

We computed the function g(u) by using Eq. (4.10), i.e. without relying on an expansion 
around the Ising FP, and by resumming the /3-functions using [3/1], [4/1], and [5/1] Pade- 
Borel approximants constrained to have a zero at v = vj = 23.56. The results up tou ^ —20 
would not be distinguishable from the quadratic approximation shown in Fig. 5. 

Let us now consider Cy.{y). This function can be computed directly by using Eqs. (4.14) 
and (4.19). They provide as a function of the variable a. In order to compute the relation 
between a and y, we need to determine the small-cr behavior of C^. We write 

Cx = l + Ecn(T", (B32) 

n=2 

and, as for g{u), we compute perturbative series Cn{v) such that c„ = Cn{vi). By resumming 
these expansions we obtain 

C2 = 0.0052(12), C3 = C3C2^^^ = -4(2). (B33) 

The variable y defined by the normalization condition (4.18) is related to o" by y = —'^2"^^ — 
-0.072(8) (7. 
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